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We have made radio- and millimeter-wavelength interferometric ob- 
servations of a variety of molecular absorption lines toward the gravi- 
ty \ tational lens systems B0218+357 and PKS 1830—211. The absorption 
^\ ■ occurs in the lensing galaxies at redshifts of 0.685 and 0.89, respectively. 

The high spatial resolution of our VLA and VLBA observations allows 
O imaging of the background continuum emission and the absorption dis- 

tribution. Our multi-transition studies yield estimates of the cosmic 
■ microwave background temperature at z = 0.89 and determinations of 

55 \ molecular and isotopic abundance ratios, which allow meaningful corn- 

ed ■ parisons with Galactic molecular clouds. 
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C3 ' 1. Introduction 



The discovery by Wiklind & Combes of high-redshift molecular absorption to- 
ward four radio sources provides a novel approach of studying the dense inter- 
stellar medium of galaxies at cosmological distances. In two of the detected 
sources, 1504+377 and 1413+135, the absorption takes place in the host galaxy 
of the active galactic nucleus (AGN) providing the background continuum emis- 
sion (Wiklind & Combes 1994, 1996a). In the other two systems, B0218+357 
and PKS 1830—211, the absorption occurs in galaxies that act as gravitational 
lenses magnifying the radio- and millimeter- wavelength emission of a background 
AGN. In both of the latter sources the absorbing material, which is, respectively, 
at redshifts near 0.685 and 0.89, has a remarkably high molecular hydrogen col- 
umn density exceeding 10 22 cm~ 2 (Wiklind Sz Combes 1995, 1996b; Menten 
& Reid 1996). Particularly toward PKS 1830—211, one observes a wide va- 
riety of absorption lines from a series of high dipole moment species that are 
well-known constituents of molecular clouds within the Milky Way (Wiklind &; 
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Combes 1996b, 1998). Here we report interferometric observations targeted at 
understanding the spatial distribution and chemical composition of the absorb- 
ing clouds observed toward B0218+357 and PKS 1830-211. 

2. Observations 

We have used the VLA to search for absorption from a number of molecules 
toward B0218+357 and PKS 1830-211. In particular, the J = 1 - rotational 
ground state lines from the high dipole moment species HCO + , HCN, and HNC, 
which are detected toward many Milky Way molecular clouds, have rest frequen- 
cies in the 3 mm band that are redshifted into the 38-49 GHz range covered by 
the VLA's "new" highest frequency receivers for z ~ 0.83 - 1.3. At the time 
of our observations, 13 VLA antennas were equipped with 38-49 GHz receivers. 
Toward PKS 1830—211, we detected absorption from the J = 1 — lines of the 
three molecules mentioned, as well as from the analogous transitions of their 13 C- 
substituted isotopomers. In addition, we detected absorption from the C3H2, 
HC3N, H2CO, and C2H molecules. Our VLA observations allow for accurate 
optical depth determinations. Moreover, since some of our observations used 
the highest resolution A configuration, we are able to image the PKS 1830—211 
system at 0''l resolution. The observations described here are part of an ex- 
tensive program of studying molecular absorption toward PKS 1830—211 and 
B0218+357 that also includes observations with the IRAM interferometer, the 
Effelsberg 100 m telescope, and the Very Long Baseline Array (VLBA). 

3. PKS 1830-211: Morphology and Spatial Extent of the Absorp- 
tion 

The upper left panel of Fig. 1 shows a 47.0 GHz VLA map of the PKS 1830-211 
gravitational lens system made with a resolution of 0''1. Two compact sources, in 
the following referred to as the NE and SW images, with weak extended "tails" 
are seen, displaying the inversion symmetry expected from gravitational lensing 
of a core/jet source. With decreasing frequency, the jet components grow larger 
(Rao & Subrahmanyan 1988; Subrahmanyan et al. 1990) and form an Einstein 
ring (Jauncey et al. 1991), which is not resolved in our 0"75 resolution 7.7 GHz 
map shown in the lower left panel of Fig. 1. 

Also shown in Fig. 1 are spectra of the HCN J = 1 — line made toward 
both images. Strong absorption with a peak optical depth ^ 2 is observed 
at zero velocity (assuming z = 0.88582) toward the compact core of the SW 
image, while much weaker absorption is observed in the same transition toward 
the NE image at a 147 km s _1 lower velocity. This general distribution of the 
HCN absorption is consistent with the results obtained for the 2—1 line of 
this molecule by Frye, Welch, & Broadhurst (1997) and Wiklind h Combes 
(1998). Our higher resolution VLA observations, together with VLBA data on 
the HC 3 N J = 5 - 4 line (redshifted to 24.1 GHz; Carilli et al. 1998) allow 
for a more detailed analysis of the spatial structure of the molecular absorption 
toward PKS 1830-211(SW). As discussed by Carilli et al. (1998), the HC 3 N 
absorption is arising from a region of size larger than 2.5 milliarcseconds or 
lO/i -1 pc, where h is the Hubble constant in units of 100 km s" 1 Mpc~ x and 
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Figure 1. VLA continuum images and molecular absorption spectra 
of PKS 1830—211. The upper and lower left panels show continuum 
maps made with the VLA at frequencies of 47 and 7.7 GHz, respec- 
tively. Note that the 7.7 GHz map shows a larger field. The 47 and 7.7 
GHz maps were restored with circular beams, indicated in the lower 
left map corners, of size 0''l and 0''75 (FWHM), respectively. The cross 
and dotted ellipse in the 47 GHz map indicate the centroid position, 
orientation, and core size of the oblate spheroidal lens mass distribution 
assumed in the lens model of Nair et al. 1993, which has a core radius 
of 0"34 and an eccentricity of 0.63. The upper two panels on the right 
hand side show absorption spectra of the HCN 1 — rotational tran- 
sition taken toward the north-eastern (NE) and south-western (SW) 
images, respectively, while the lower right panel shows a spectrum of 
the 2n — 2i2 line of H2CO toward the SW image. The velocity scales 
of all spectra are in the heliocentric system with zero velocity corre- 
sponding to z = 0.88582. 
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Figure 2. Spectra of the redshifted J = 1 — transition from 
the main isotopomers of HCN, HNC, and HCO + observed toward 
PKS 1830—211 (SW) are shown in the upper row, while spectra of 
the same transition from the 13 C substituted isotopomers are shown in 
the lower row. The velocity scales of all spectra are in the heliocentric 
system with zero velocity corresponding to z = 0.88582. 



qo = 0.5 is assumed. An upper limit of ~ 600/i pc on the size of the absorbing 
cloud follows from the fact that we do not detect HCN absorption toward the 
tail of SW image with an upper limit on the optical depth of 0.3. 

As discussed by Wiklind & Combes (1998), the observed difference in ab- 
sorption velocities between the SW and NE images allows for tight constraints 
on the lens mass distribution and, thus, lensing models of PKS 1830—211 (see 
Kochanek & Narayan 1992; Nair, Narasimha, & Rao 1993). Wiklind & Combes 
(1998) derive Vq 220\^D km s _1 for the rotation velocity, Vq, of the z ~ 
0.89 lensing galaxy. The "effective distance" D, measured in Gpc, is equal to 
D\D s j Dig] where Di is the distance between the observer and the lens, D s the 
distance between the observer and the source, and D[ s the distance between 
the lens and the source. All distances are angular diameter distances. Very re- 
cently, a value of 2.507 has been measured for the redshift of the lensed source in 
the PKS 1830—211 system (see the report by Lidman 1998), yielding V$ ps 366 
km s _1 and an inclination angle, i, of 16° (Hq = 75 km s" 1 Mpc -1 and qo = 0.5 
are assumed throughout this paper). This indicates that the lensing galaxy is a 
massive early-type spiral. There is, however, the possible complication of mul- 
tiple lensing, which might be expected, given the presence of an HI absorption 
system at a redshift of 0.1927 (Lovell et al. 1996). 



4. Astrochemistry at z = 0.88582 

4.1. Gas Excitation and Temperature of the Cosmic Microwave Back- 
ground 

Figs. 1, 2, and 3 show examples of absorption line spectra toward the PKS 
1830—211 (SW) source resulting from our VLA observations. For some species 
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Table 1. 



Molecules toward PKS1830-211 and Galactic Clouds 



Species 


H 


N(PKS 1830-211) 


AT(X)/7V(CO)xl0 6 












(D) 


(cm- 2 ) 


PKS 1830-211 TMC-1 


L 183 


3C111 


NRAO 150 


CO 


0.11 


3.0xl0 18 


1 


1 


1 


1 


1 


cs 


1.96 


4.8xl0 14 


160 


130 


9 






HCO+ 


3.30 


4.0xl0 14 


130 


110 


100 


130 


180-1400 


HCN 


2.98 


3.2xl0 14 


110 


130 


40 




400-2000 


HNC 


2.70 


1.2xl0 14 


40 


250 


<50 






N 2 H+ 


3.40 


5.4xl0 13 


18 


10 


~10 






CCH 


0.8 


7.9xl0 14 


260 


>100 


<100 






H 2 CO 


2.33 


3.8xl0 14 


130 


250 


250 


220 


590-1200 


C 3 H 2 


3.4 


6.7xl0 12 


2 


130 


3 


56 


190 


HC3N 


3.72 


8.1xl0 12 


3 


75 


1 






OH 


1.67 


4.1xl0 15 


1400 


3800 


940 


250 





The first, second, and third columns list, respectively, molecular species detected toward 
PKS 1830—211, their dipole moments (in Debyes), and column densities in the z — 0.88582 
absorbing cloud measured toward PKS 1830— 211(SW). Column densities of species listed in 
italics are not derived from our data but are taken from Wiklind & Combes 1996b and, for 
CO and OH, from Gerin et al. 1997 and Chengalur, de Bruyn, & Narasimha (this volume), 
respectively. The remaining columns list for PKS 1830—211 and four well-studied Galactic 
molecular clouds the column densities measured for the various species multiplied by 10 6 and 
divided by the column density of CO in the clouds in question. The Galactic clouds considered 
are the dark clouds TMC-1 and L 183 (also known as L 134 N; see Irvine et al. 1987 and 
van Dishoeck et al. 1993 and references therein for abundances) and the diffuse clouds seen in 
absorption against the extragalactic radio sources 3C111 and NRAO 150 (Liszt & Lucas 1995, 
1998; Lucas & Liszt 1996; and Cox, Glisten, & Henkel 1988). For NRAO 150 the range of 
values observed for the various absorption components occurring at different radial velocities 
is given, except for the case of C3H2. 
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Figure 3. Redshifted absorption spectra of various hyperfine compo- 
nents of the N = 1 — transition of CCH (left hand panels), the HC3N 
J = 3 — 2 and 5 — 4 transitions (middle panels), and the 2^ — loi 
transition of C3H2 (right hand panel) observed with the VLA toward 
PKS 1830—211. The velocity scales of all spectra are in the heliocentric 
system with zero velocity corresponding to 2 = 0.88582. 



we can compare our measured optical depths with the values determined by 
Wiklind & Combes (1996b) for higher excitation lines of the same species. We 
find that to within the uncertainties of a few K, the thus determined rotation 
temperatures, T rot , are consistent with the value of the cosmic microwave back- 
ground temperature, T cm b, at z = 0.89. All of the species observed by us have 
high dipole moments (see Table 1). Even for the relatively low frequency HC3N 
lines, the critical densities necessary for thermalization of the level populations 
at or near the kinetic temperature, Tk in , are ^ 10 4 cm -3 . Therefore, the low ro- 
tation temperatures indicate subthermal excitation due to a relatively low total 
density of the absorbing cloud rather than low kinetic temperatures, which for 
molecular clouds within the Milky Way are generally greater than 10 K. 

Given that T rot = T cm ^, the accurate optical depth determinations afforded 
by our interferometer measurements allow for a meaningful estimate of the cos- 
mic microwave background temperature at z = 0.89, for which big bang theory 
predicts a value of (1+z) 2.73 K = 5.14 K. In particular, from the HC3N J = 3—2 
and 5 — 4 spectra shown in Fig. 3 we derive T rot = 4.51J' g K. We note that the 
errors quoted for the HC3N rotation temperature are formal uncertainties and 
do not take into consideration systematic effects such as variations in the source 
covering factor between 14.5 and 24.1 GHz, the frequencies of the redshifted 
3 — 2 and 5 — 4 lines. 

4.2. Molecular Abundances — Comparison with Galactic Clouds 

From our measured optical depths we calculate the total column densities for 
HCO+, HCN, HNC, CCH, H 2 CO, C 3 H 2 , and HC 3 N listed in Table 1, assuming 
that the level populations are thermalized at T cm h — 5.14 K. Table 1 also shows 
a comparison of the abundances of molecules detected toward PKS 1830—211 
with the values measured toward four Galactic molecular clouds. Listed are the 
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abundances of the various species (multiplied by 10 6 ) relative to the CO column 
density. The Galactic clouds in question are the cold, dense dark clouds TMC-1 
and L 183 and the more diffuse clouds seen against the extragalactic radio sources 
3C111 and NRAO 150. The dark clouds have CO column densities, JV(CO), of 
order a few times 10 18 cm" 2 and [CO/H2] abundance ratios «8x 1CT 5 (Irvine 
et al. 1985; van Dishoeck et al. 1993 and references therein); while Liszt & Lu- 
cas (1998) determine iV(CO) = 9 x 10 16 cm" 2 for 3C111 and iV(CO) values in 
the range (0.25 - 6.6) xlO 15 cm" 2 for the various radial velocity components of 
the multiple-cloud absorption system observed toward NRAO 150. The 3C111 
absorber and most of the absorbing clouds seen toward NRAO 150 are quite 
dense compared to "typical" diffuse clouds studied, e.g., by ultraviolet absorp- 
tion spectroscopy and probably have [CO/H2] abundance ratios exceeding 10" 5 . 
By comparison, for the diffuse cloud sample discussed by Federman et al. (1994; 
see also Liszt & Lucas 1998) [CO/H2] abundance ratios range from a few times 
10" 8 to a few times 10" 6 . 

It is clear from Table 1 that the PKS 1830—211 absorbing cloud has a total 
column density and molecular abundances that are quite similar to the values 
found in Galactic dark clouds (note that C3H2 and HC3N are anomalously over- 
abundant in TMC-1 by Galactic molecular cloud standards). The z = 0.88582 
cloud does not show the surprisingly high abundances of HCO + , HCN, and 
H2CO relative to CO found in the more diffuse Galactic clouds seen in absorption 
toward 3C111 and NRAO 150 and other sources, which are, particularly in the 
case of H2CO, a challenge to current models of diffuse cloud chemistry (Turner 
1994; Liszt & Lucas 1995). 

While there are chemical similarities, a comparison of cloud structure in 
space and velocity reveals remarkable differences between the z = 0.88582 ab- 
sorbing cloud and typical Galactic dark clouds such as TMC-1 or L 183. As- 
sume that for the PKS 1830—211 (SW) cloud the abundance ratios of, say, the 
widespread CO, HCO + , and HCN molecules relative to H2 are similar to the 
values found for Galactic dark clouds. This implies a total molecular hydrogen 
column density of a few times 10 22 cm -2 . Assuming further that the molecular 
hydrogen density is smaller than 10 4 cm" 3 , as implied by the subthermal excita- 
tion of all of the high dipole moment species, we derive a pathlength larger than 
~ 1 pc along the line of sight, which is to be compared with the 13 pc lower limit 
of the source size in the plane of the sky indicated by the VLBA results of Car- 
illi et al. (1998) discussed above. To compare with Galactic molecular clouds, 
it is instructive to take as a reference source the nearby and well-studied Tau- 
rus/Auriga dark cloud complex (see, e.g., Cernicharo 1991), which has a total 
extent of ~ 50 pc, while individual clouds within that complex, such as TMC-1, 
have sizes of one to a few pc. The line widths in such clouds are < 1 km s" 1 , 
and radial velocities are within ~ 1 km s" 1 of a median value, with only a few 
km s" 1 variation over the whole extent of the complex. In contrast, from, e.g., 
the high quality spectrum of the optically thin H2CO absorption line (see Fig. 
1), we find that the z = 0.88582 absorption covers a total velocity range of ^ 35 
km s" 1 (FWZP), which is significantly larger than the velocity spread even of 
a Galactic giant molecular cloud (GMC), with the exception of the remarkable 
GMCs found within ~ 200 pc of the Galactic center. The latter clouds have ve- 
locity widths that are comparable to that of the z = 0.88582 absorber. We note 
that compared to the distance of these GMCs to the Galactic center, the cloud 
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absorbing PKS 1830—211 (SW) is located much further away from the center 
of the lensing galaxy, 1.8 kpc in the context of the model of Nair et al. (1993) 
and 3 kpc in the model by Kochanek & Narayan (1992). This possibly indicates 
significant differences in molecular cloud structure and composition between the 
z = 0.88582 lensing galaxy of PKS 1830-211 and the Milky Way. 

4.3. Isotopic Abundance Ratios 

Fig. 2 shows spectra of the redshifted J = 1 — transitions of HCN, HNC, 
and HCO + together with spectra of the analogous transitions from the 13 C- 
substituted isotopomers. Of the three main isotopic lines, the HNC line has 
the lowest optical depth (ps 1.2), affording, together with the optically thin 
H 13 CN line, a meaningful estimate of the [ 12 C/ 13 C] ratio; we find a value of 
ps 35. The [ 12 C/ 13 C] ratios determined from HCN and HCO + are consistent 
with this number. While the [ 12 C/ 13 C] ratio in the z = 0.88582 cloud is smaller 
than the values found in the solar system and the local interstellar medium, 
[ 12 C/ 13 C] ratios around 50 are measured for inner Galaxy clouds within 4 kpc 
of the Galactic center, while the Galactic center clouds themselves have even 
smaller values around 20 (see Wilson <fe Matteucci 1992; Wilson & Rood 1994). 

13 C is only produced in low and intermediate mass stars, while 12 C is also 
produced in massive stars. Thus, the [ 12 C/ 13 C] ratio is expected to decrease 
in time and with increasing stellar processing (see the discussion of Wilson h, 
Matteucci 1992). The relatively low value measured at z = 0.88582 may thus 
be of interest in the context of chemical evolution studies. However, we note 
that we only sample one line of sight in the lensing galaxy, which might have an 
atypical chemical history. 

We also searched unsuccessfully for absorption the DCN (1 — 0) transi- 
tion. Our upper limit and other searches for deuterated molecules toward 
PKS 1830—211 are discussed by Shah et al. (this volume). 

5. A Multi- Transition Study of Formaldehyde at z = 0.68466 toward 
B0218+357 

Formaldehyde (H2CO) was first observed toward the Einstein ring B0218+357 
in its 2 cm rest wavelength 2n — 2±2 If -doublet transition by Menten & Reid 
(1996), who found absorption toward the south-western of the two compact 
images of the background source, labeled A in Fig. 4 (see also Patnaik et al. 
1993 and Patnaik & Porcas, this volume, for a description of the radio continuum 
morphology). The H2CO absorption occurs at z = 0.68466, the same redshift 
at which Wiklind & Combes (1995) detect absorption from other molecules. 

The populations of the lio — In and 2\\ — 2\i if -doublet transitions at 6 and 
2 cm wavelength, respectively, are sensitively dependent on the H2 density and 
these lines are usually anti-inverted. Statistical equilibrium calculations yield ex- 
citation temperatures and optical depths, so that comparisons with observations 
in principle allow meaningful estimates of the molecular hydrogen density (see, 
e.g., Henkel, Walmsley, & Wilson 1980). While both the 6 and the 2 cm transi- 
tions have now been detected toward B0218+357 (Fig. 4), their interpretation 
is complicated by the fact that the extent and intensity of the background con- 
tinuum emission is strongly dependent on the observing frequency (see Patnaik 



8 



S,ine(Jy) 

o 

-0.003 
-0.006 



-0.002 
-0.004 



-0.04 
-0.08 

-0.12 



-0.05 

-0.1 
-0.15 



H 2 CO toward B021 8 + 357 




-30 -20 -10 10 20 30 
Velocity Offset (km s 1 ) 



B02 18 + 357 ' 8.6, GHz- 




_i i i_ 



2 10-1 

East Offset (arcsec) 



H 2 C0 Energy Level Diagram 
40 r 



35 -J K 
30 
25 
20 
15 
10 

5 





i i . 



& 12 

150.5 GHz 

l,n_ 
Ul 



1, 



140.8 GHz 
= — 4.8 GHz 



K Q = 
Para 



K a = 1 
Ortho 



Figure 4. The left four panels show redshifted absorption spectra 
of several H2CO transitions toward B0218+357. The top and bottom 
spectra were taken with the Effelsberg 100 m telescope and the IRAM 
interferometer, respectively, while the other two spectra were observed 
with the VLA. The labels in the left bottom corner of each spectrum 
list the transition and rest frequency. The O or P indicates whether the 
line in question is from ortho- or para-H^CO. The velocity scales of all 
spectra are in the heliocentric system with zero velocity corresponding 
to z = 0.68466. The upper right panel shows a VLA continuum map of 
the B0218+357 gravitational lens system (from Menten & Reid 1996) 
made at 8.6 GHz, the frequency of the redshifted 2n — 2i2 transition. 
The H2CO absorption occurs to the southwestern compact source A. 
The lower right panel shows an excerpt of the H2CO rotational energy 
level diagram. Transitions that we measured toward B0218+357 are 
shown as arrows, with rest frequencies indicated. 
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& Porcas, this volume). Modeling of the H 2 CO excitation will be greatly aided 
by VLBI observations and determinations of the excitation temperatures of the 
6 cm and 2 cm lines, which can be extracted from observations of millimeter- 
wavelength lines interconnecting the energy levels of the centimeter lines (see Fig. 
4). For this reason, we have also observed a number of millimeter- wavelength 
H2CO transitions toward H2CO. While spectra for several lines are shown in Fig. 
4, the results of our model calculations will be presented in a future publication. 

Since we have observed both ortho- and para-H^CO lines, our observations 
allow an estimate of the ortho-to-para ratio, which should be 3 if the H2CO was 
formed in warm gas (Tki n 3> 15 K) or on hot dust grain surfaces and ps 1 for 
formation in cold gas or on cold grain surfaces. 

From a preliminary analysis of the loi — Ooo para- and the 2i2 — In ortho- 
transition (see Fig. 4), we derive an ortho-to-para H2CO ratio of 3.1 ± 0.2. 
Systematic uncertainties of this value are likely to be small, since, contrary to the 
case of the low frequency if-doublet transitions mentioned above, the covering 
factors of the two millimeter-wavelength transitions used in its derivation should 
be very similar. 

Acknowledgments. The National Radio Astronomy Observatory is a fa- 
cility of the National Science Foundation operated under cooperative agreement 
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